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We use sampling techniques to find robust constraints on the masses of a possible fourth sequential 
fermion generation from electroweak oblique variables. We find that in the case of a light (115 GeV) 
Higgs from a single electroweak symmetry breaking doublet, inverted mass hierarchies are possible 
for both quarks and leptons, but a mass splitting more than Mw in the quark sector is unlikely. We 
also find constraints in the case of a heavy (600 GeV) Higgs in a single doublet model. As recent 
data from the Large Hadron Collider hints at the existence of a resonance at 124.5 GeV and a single 
Higgs doublet at that mass is inconsistent with a fourth fermion generation, we examine a type II 
two Higgs doublet model. In this model, there are ranges of parameter space where the Higgs sector 
can potentially counteract the effects of the fourth generation. Even so, we find that such scenarios 
produce qualitatively similar fermion mass distributions. 

PACS numbers: 14.65.Jk 



I. INTRODUCTION 

Adding a sequential fourth generation of fermions (4G) 
is one of the simplest possible extensions to the Standard 
Model. Indeed, although the width of the Z limits the 
number of active light neutrinos to three, there can in 
principle be a fourth neutrino generation which is much 
heavier. A recent and extensive literature examines the 
impact of how 4G would reduce tensions in recent mea- 
surements in the b sector and create distinctive phenom- 
ena in kaon decays [l|, , as well as providing a potential 
scenario for the baryon asymmetry of the universe [3( . As 
a more top down motivation, simple string constructions 
often lead to toy models with an even number of gen- 
erations. Of course, achieving three chiral generations 
is also possible in a wide class of examples, and some 
stringy models of flavor physics predict that more than 
three generations would be inconsistent with the mea- 
sured three generation quark mixing matrix Q. More 
generally, one can view the 4G scenario as a simple tem- 
plate for scenarios of physics beyond the Standard Model 
in which states of some extra sector receive a mass pro- 
portional to the Higgs vev. 

In light of these considerations, it is clearly of inter- 
est to study the viability of the 4G scenario. In addi- 
tion to the possibility of direct detection of such states, 
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the contributions of these additional states enter as loop 
corrections to various Standard Model processes. For 
example, a fourth generation tends to produce a positive 
contribution to both S and T, the oblique electroweak pa- 
rameters. By contrast, in a single Higgs doublet model, 
increasing the mass of the Higgs generates a positive con- 
tribution to S and a negative contribution to T. Thus, 
while cancellation for the T parameter is possible, the 
contributions to the S parameter typically move in the 
same (positive) direction, though mass hierarchies in the 
fourth generation can reduce the size of this contribution. 
The extra generation also affects the phenomenology of 
the Higgs, leading to an increase in T(h — > gg), and a 
decrease in T(h — > 77). 

Though less well-studied, even simple extensions of the 
Higgs sector can counteract (or exacerbate) some of the 
effects of a chiral fourth generation with an appropriate 
tuning of parameters. For example, in two Higgs doublet 
models (2HDM), the contributions to S and T can have 
either sign (see e.g. 0-0] )• General values of the Higgs 
mixing angles also allow for changes in T(h — > gg) and 
T(h — > 77), independently, relative to the 4G scenario. 

In this paper we study the available parameter space 
for the 4G scenario and its extension to 2HDM models. 
The full parameter space of 4G is too large for easy visu- 
alization, but much of what we need to know to under- 
stand existing experimental results and to inform future 
searches can be expressed with two pairs of numbers: 
the two quark masses and the two lepton masses. We 
therefore seek, by sampling this four-parameter space and 
comparing the samples with constraints on electroweak 
oblique parameters to determine the most likely mass 
spectrum for 4G, should it exist. Similar earlier analyses 
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of this type may be found in [il-fioj , but new experimen- 
tal data has appeared since these publications. Other 
similar studies have appeared recently fliT-[l3j]. 

Exclusion limits on the mass of a Standard Model- 
like Higgs impose additional constraints on the 4G sce- 
nario. In 4G, the gluon fusion production cross sec- 
tion for the Higgs is markedly increased over the three 
Both LHC collaborations 
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generation scenario, 
dependently exclude, using a combination of channels, 
the range 120 GeV < M h < 600 GeV when there is a 



fourth generation. The LEP II lower limit of 114.4 GeV 
is independent of the number of fermion generations [15l | . 
Because a fourth generation of fermions contributes to 
T roughly quadratically with m{ti) and 771(64), and be- 
cause a large T corresponds to a large m(h), values of 
rn(h) as hig h as 1 TeV are allowed by electroweak con- 
straints jlOl [l6| in 4G. However, studies of the stabil- 
ity and triviality bounds on m(h) in 4G [l7| prohibit 
m(h) > 700 GeV unless there is also some other new 
phenomena on a scale below 2 TeV. 



3 



150- 
100 ^ 

50 E- 
o E- 

-50 E- 
-100 E- 
-150 E- 
-200^ 




0.016 
= 0.014 
0.012 
0.010 
0.008 
0.006 
^0.004 
Ho.002 



100 150 200 
m(T 4 ) - m(v 4 ) 



■o.ooo 



200 
150 
100 

50 


-50 
-100 
-150 



— I I I I I I I — 











4 



0.025 
0.020 
0.015 
0.010 
0.005 



100 150 200 
m(x 4 ) - m(v 4 ) 



FIG. 3: Contour plots of the probability densities for quark vs. lepton mass splitting in the 4G scenario. Left: with m(h) = 
115 GeV; Right: with m(h) = 600 GeV. The boxes mark the areas where the magnitude of the mass splittings is less than Mw- 
All scales are in GeV; probability densities have been normalized and each bin is 10 GeV x 10 GeV. 



Most recently, there are "hints" of a Higgs with 
mass [3 of 124.5 ± 0.8 GeV from the LHC [H and sup- 
porting evidence from the Tevatron [20j. The hint is 
strongest in the channel gg — > h — > 77, where the AT- 
LAS experiment reports an excess above background of 
2.8 standard deviations. The statistical significance of 
these results is not enough to declare discovery or even 
strong evidence for a Higgs, but is strong enough to pro- 
voke discussion. This mass is within the bounds ruled 
out by the LHC when supposing a fourth fermion gener- 
ation. To leading 0(GpnVf), it is possible to retain 4G if 
one supposes only the 77 channel's hint remains signif- 
icant with the addition of more data, but including ex- 
act next-to-leading order electroweak corrections makes 
this difficult [ll], [22| . Consequently, the 4G hypothesis is 
valid only if (a) the hints turn out to be statistical fluc- 
tuations or (b) the hints are due to something beyond a 
single Higgs doublet, such as a two Higgs doublet model. 

The rest of this paper is organized as follows. In Sec- 
tion|Hjwe treat the 4G single Higgs doublet case. We give 
results for m(h) = 115 GeV (the "baseline" scenario) 
and for m(h) = 600 GeV (the "high-mass" scenario). 
The baseline scenario is appropriate for the case con- 
sidered in [U [22[ ; the difference in our results between 
m(h) = 115 GeV and 124.5 GeV is small. In Section |HT] 
we extend the analysis to consider a Type II model with 
its parameters adjusted to match the 124.5 GeV hint. 
Section ITVl provides a summary. 

II. SINGLE HIGGS DOUBLET SCENARIOS 
A. Method 

We have updated constraints on the oblique elec- 
troweak [HI parameters S, T and U, as found by the 
Global Analysis of Particle Properties (GAPP) [IJj using 
data available in October 2011. In our sampling proce- 
dure, each sample is assigned a weight corresponding to 



the probability density function p — p(S, T, U) for these 
three parameters. We employ the one-loop contributions 
to the oblique parameters, assuming small mixing with 
the extra family, as in Q . See [25| for some recent discus- 
sion of the more general case of potentially large mixing 
effects. 

The sampling distribution in this type of analysis plays 
the role of a Bayesian prior; we are taking the probability 
p of a specific value for S, T and U given an assumed set 
of four fermion masses, and weighting it in our result as 
the probability density created by our sampling of the 
fermion spectrum. We interpret the result as a probabil- 
ity density function for the fermion mass spectrum, but 
that interpretation is only valid in the context of that as- 
sumed sampling distribution. The peril in this process - 
the validity of the assumed prior - thus has the advantage 
of requiring explicit description. 

We draw 50 million uniformly distributed samples in 
the fermion mass spectrum with lower bounds set by 
direct experimental constraints described below. The 
upper bound is limited by unitarity arguments [26[ to 
500 GeV, but this is a rough bound and we raise it to 
700 GeV for clarity in the resulting figures. 

The lower bound on the sampled 771(1/4) mass range is, 
in our baseline scenario, 777(7/4) = 90.3 GeV from LEP 
II [27J. This limit is the weakest of the limits obtained 
under the assumption of 7/4 decay to each of the three 
known charged leptons; if 777(7/4) > 777^4), then we would 
obtain a stronger limit. The lower bound on the sampled 
777(74) range is 100.8 GeV; again, this is the weakest limit 
obtained in all the possible decay scenarios. These re- 
sults are therefore robust against all assumptions about 
the lepton mass hierarchies. On the other hand, lep- 
ton mixing parameters are important considerations in 
searches for the leptons of 4G at the LHC which have 
been discussed [38| but have not yet been carried out. 

Obtaining robust lower bounds on 4G quark masses 
and mixing angles is a little more complicated. Dramatic 
results [58[ from the LHC are indeed available [28[, and 
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new ones are appearing constantly. The CMS collabora- 
tion lias searched for: 



• 6464 tW-tW+ -> (bW+W~)(bW-W+) 
same-ch arg e leptons and trileptons in a 4.6 fb 
sample [29[ , obtaining a limit of 600 GeV. 



with 

-1 



• both i 4 and 64 using a simplified model with a 
range of final states, all containing 2 6 quarks, in 
f.f ftT 1 of data All of the diagrams consid- 
ered have 64 — > tW or £4 — > bW. Lower limits of 
480 - 540 GeV were obtained. 

• pair produced t± in the "lepton with jets" channel, 
wherein a decay to bW having the same signature as 
a it event but with a different primary quark mass 
is sought. The analysis reconstructed m{ti) in each 
event. A 560 GeV lower limit was found using only 
4.6 ftT 1 of data [H|. 

• pair produced tn in the "dilepton" channel, wherein 
also a decay to bW having a top-quark signature 
but different mass is sought. A weaker constraint 
than that which was obtained in the "lepton with 



jets" analysis, 422 GeV, was found using f .f fb 
data 



of 
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The ATLAS collaboration has searched for: 

• pair produced 64 — > tW in 34 pb -1 of data [3J|, as 
part of an inclusive search for exotic production of 
the same-charge dilepton signature. 

• pair produced £4 or 64 decaying to Wq, where 
q = u,d,s, or b, appearing with opposite-charge 
dilcptons and missing transverse momentum in 
37 pb -1 of data |34j |. An approximate event recon- 
struction is done. The resulting limit is 777 (£4) = 
m(6 4 ) > 270 GeV. 

• 6464" -> tW~tW+ -> (bW + W-){bW~W+) with 
one lepton, at least six jets, and large missing mo- 
mentum transverse to the beamline on a 1.0 fb -1 
sample, obtaining |35J a limit of 480 GeV. 

• pair produced £4 or 64 appearing with same-charge 
dilcptons, large missing transverse momentum, and 
at least two jets in 1.0 fb -1 of data [3g|. A limit of 
771(64) > 450 GeV was obtained. 

See [37| for recent searches of more exotic fermions. 

These search results, while impressive, are all built 
upon specific decay, i.e. CKM mixing angle, assumptions. 
With the exception of 34|, mixing of the fourth genera- 
tion into anything other than the third generation is not 
considered. Furthermore, t± — > b^W* (or in an inverted 
hierarchy, 64 — > t^W*) will be an additional contribution 
to 64 (or /j 4 ) production which will not necessarily appear 
in any specific signature as a result of the W* products; 
the contribution from this channel can be significant if 
the mass splitting is small. 



Additionally, there are constraints on the possible mix- 
ing parameters. For example, the mixing parameters for 
the quark sector may be constrained [2|, |39| with data 
from neutral mesons, the 6 — > 57 transition, existing 
constraints on the three-generation quark mixing matrix 
and limits on Br(£>s ~~ )• Reference 3!) concludes 
that large mixings of the fourth generation with the three 
known generations are not ruled out, but [4Ct |41|. which 
considers constraints from corrections to the Z — >bb ver- 
tex from a fourth generation conclude that these mixings 
could be comparable to Cabibbo mixing. The quark mix- 
ing matrix can also be constrained with precision elec- 
troweak data and D° — D° mixing [42j . In any case how- 
ever, there is the possibility that 4G fermions could decay 
to either third or lower generation fermions with varying 
branching ratios. 

A method for producing experimental limits that are 
mixing-angle independent [59j and that allows for the 
contributions of both 4G quarks to any particular sig- 
nature was applied to the results of CDF searches [43j j . 
resulting in lower limits of ~ 280 GeV for 777(64) and 
~ 290 GeV for m{ri). We use these lower but mixing 
independent values here while strongly advocating the 
application of these techniques to the more recent LHC 
results. Such an analysis could soon sharply constrain or 
even rule out the 4G hypothesis. 



B. Results 

Figures [T] and show the lepton and quark mass spec- 
tra in our baseline and high m(h) scenarios. In these 
and similar Figures, the color for each bin represents 
a probability density integrated over the bin, and nor- 
malized so as to give unit probability when summed 
over the entire plot. For the baseline (high mass) case, 
\m(ti) - 777(64)1 < M w in over 99% (90%) of our sam- 
ples; transitions between 4G quarks will produce off-shell 
W bosons. The lepton mass splitting is less than M\y 
with probability 69% (24%). Normal mass hierarchies 
are more likely than not, but by no means certain; in the 
lepton sector the probability of a normal mass hierarchy 
is 70% (93%), and in the quark sector, it is 59% (69%). 

Masses just over the existing limit for the leptons are 
heavily favored, and this tendency is greater in the high 
m(h) scenario. Being able to predict this parameter rel- 
atively precisely makes it a valuable target for future 
searches. 

In Figure [3] we show the lepton and quark mass split- 
tings. We see that, with perhaps a two-fold ambiguity, 
the mass splittings in the two sectors are tightly related. 

Carpenter and Rajaraman [44( revisited the LEP II re- 
sults in a scenario with both left- and right-handed neu- 
trinos. They conclude that 777(1/4) as low as 62.1 GeV is 
possible. Some recent studies [ll], EH H2| also consider 
low values of 777(1/4). We find that lowering the bound 
on 777(7/4) to 62.1 GeV does not produce much change rel- 
ative to our baseline scenario. Figure U shows distribu- 
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tions that have the same probabilities of mass splittings 
less than M\y and the same probabilities of normal mass 
hierarchies as our baseline scenario to within about 2%. 



III. TWO HIGGS DOUBLET SCENARIO 
A. Method 



Should the hints of a Higgs boson with m(h) = 124.5 
GeV solidify with more data, the 4G hypothesis is only 
tenable if an extended electroweak symmetry breaking 
sector exists. As an example of such an extension we 
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FIG. 6: Left: Contour plots of the probability densities in 2HD4G with the mass of the lightest CP— even state m(h) = 124.5 
GeV. Quark mass splitting vs. lepton mass splitting. The box marks the area where the magnitude of the mass splittings is 
less than Mw- All scales are in GeV; probability densities have been normalized and each bin is 10 GeV x 10 GeV. Right: The 
probability density function for log 10 (tan/3) in 2HD4G with the mass of the lightest CP— even state m(h) = 124.5 GeV. 




consider a second Higgs doublet [45j in conjunction with 
a fourth sequential fermion generation (2HD4G). Two 
identical complex scalar SU(2)l doublet fields $i and 
$2, both of hypercharge Y = 1 are postulated. To for- 
bid flavor changing neutral currents, we select the Type 
II Yukawa coupling pattern, in which Q = 2/3 quarks 
couple to one doublet and Q = —1/3 quarks and 
charged leptons to the other. This restriction permits 
a Z2 symmetry to distinguish $1 from $2- We restrict 
consideration to the gauge invariant, renormalizable and 



CP— conserving potential 

V = ml^l^ + m| 2 $ 2 $ 2 + y (O^i) 2 + y (*^ 2 ) 2 
+ A 3 ($l$i)($t$ 2 ) + A 4 ($ t 1 $ 2 )(^$i) 

+ y((*i*2) 2 + (<^i) 2 ) (1) 

where all the parameters ma and Aj are real. This sys- 
tem and its vacua preserve an additional I2 symme- 
try. There are two CP— even neutral bosons, h and H 
(m(h) < m(H)) a CP— odd neutral boson, A, and the 
charged bosons H in this model. 

This model is different from the similarly-named 
"4G2HDM" model of [Hj]; however, [47| analyzed a simi- 
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lar model prior to the appearance of the 124.5 GeV hint. 
The presence of large fourth generation Yukawas can lead 
to large radiative corrections which can potentially desta- 
bilize the form of the Higgs potential. Here we assume 
that the 2HDM effective potential is stabilized by some 
effect near the TeV scale, so that we can focus on the 
resulting effective theory below the TeV scale. 

Though it is beyond the scope of this study, the com- 
bination of two Higgs doublets with a fourth sequen- 
tial fermion generation creates a rich phenomenology for 
which constraints from the kaon and b sector could be de- 
rived. For example, the coupling constants Z —> bb vertex 
will obtain corrections which depend on T4 4 &, m{t/Cj, and 
(depending on chirality) either tan/3 or its inverse; these 
contributions can be constrained experimentally. 

Two important parameters of this model are tan /3, 
the ratio of the vacuum expectation values of the two 
doublets and a, the angle which diagonalizes the mass- 
squared matrix of the CP— even bosons. Values of tan /3 
less than 1 are disfavored experimentally assuming three 
fermion generations; more generally, tan j3 > 0.3 results 
from the requirement that the top quark Yukawa cou- 
pling not exceed the perturbative limit (45| . Requiring 
perturbativity of the fourth generation Yukawa interac- 
tions can impose additional constraints. For the sake 
of generality, however, we do not impose this additional 
restriction in our scans. We sample tan/3 in a scale- 
independent way, i.e., the distribution of log(tan/3) is 
uniform. The angle a is scanned uniformly but samples 
are weighted according to the value of a as described be- 
low; the masses of the 4G fermions and the bosons H, A 
and are selected with an initially uniform distribu- 
tion. The mass of the lightest CP— even boson is set to 
124.5 GeV. For further discussion on the phenomenology 
of two Higgs doublets with a fourth fermion generation, 
including the case where v 4 is stable and contributes in- 
visible decays to either h or A, see [H]]. 

The 124.5 GeV hint is strongest in the channel gg — > 
h — > 77, where the ATLAS experiment reports an excess 
above background of 2.8 standard deviations. The second 
most significant hints are in the channels gg — > h — > VV*, 
where the ATLAS results have a significance of 2.1 and 
1.4 standard deviations for V = W, Z respectively. The 
combination of ATLAS and CMS data correspond to a 
77 production rate about 1.4 ± 0.7 times the predic- 
tion of the Standard Model [II]; for VV* , it is about 
0.8 //times the Standard model rate. 

For each scanned value of a, we calculate 

a(gg — > M r ( ft ^"V/j) anc ; a(gg — > h) ^ h ^ v Y. \ 

^ ' 1 2HD4G(tot) ' 1 2H D4G (tot) 

for the 2HD4G scenario and form a \ 2 of these values 
against these experimental values. We weight each sam- 
ple according to that x 2 - We do not consider constraints 
from decays of the H , A and H ± which are very param- 
eter dependent in the 2HD4G scenario. 

The dominant production mechanism at both the LHC 
and the Tevatron is gluon fusion through loop diagrams 
involving the colored states. In a 2HDM, this includes 



the contributions from both the t, b as well as and 
64. The Standard Model normalized cross section a = 
o'2HD4g/ ij sm of the gluon fusion production cross sec- 
tions is: 

^ _ IgCAyag) + A 1/2 (U)) - %(A 1/2 (b) + A 1/2 (b 4 ))\ 2 

\A 1/2 (tW 

(2) 

where c a and s a denote cos a and sin a respectively, and 
A-1/2 (A) is the threshold correction of a spin 1/2 particle 
X to the h — > cm vertex for a 124.5 GeV Higgs, with 
notation as in 50]. A similar expression holds for the 
Standard Model normalized decay rate h — ¥ 77. In a 
2HDM, this will include terms from loops containing W , 
t, b and r and charged fourth generation fermions, as 
well as a contribution from H ± , which all depend on the 
mixing angles. The total width of the Higgs in 2HD4G 
including the mixing angle dependence is fixed by similar 
considerations. Much as in [511 ] . the overall normalization 
can be extracted from the recently updated values for the 
Standard Model 124.5 GeV Higgs partial widths [E3] by 
including the mixing angle dependence and contribution 
from extra states in the various 2HD4G partial widths. 

Constraints on two doublet models are readily avail- 
able [HH through the package 2HDMC. We observe 
the constraints of tree- level unitarity [55| , perturbativity 
(i.e., the magnitudes of all the quartic Higgs couplings 
must be less than 47r), and the absence of runaway direc- 
tions, as implemented in 2HDMC. Contributions to the 
oblique electroweak parameters [57J are also provided as 
part of 2HDMC. 



B. Results 

Figures [5] and [5] show the lepton and quark mass spec- 
tra in our two Higgs doublet scenario. The quark (lepton) 
mass splittings are less than M\y in 99% (65%) of our 
samples; normal mass hierarchies in the quark (lepton) 
sector occur with a probability of 59% (72%). 

Low values of tan /3 are likely in 2HD4G; in Figure |H1 
tan /3 < 1 in 46% of the final probability density func- 
tion. Figure [7] shows the distribution of Higgs boson 
masses. There is a strong correlation between the masses 
m(P ± ) and m(A) largely but not entirely created by re- 
quiring VV* as well as 77 production to be in agreement 
with experiment. It is amusing to note that the most 
likely values for the mass of the second CP— even boson 
are just over 124.5 GeV, and masses corresponding to a 
small excess in the 4£ channel at 240 GeV are not im- 
probable. 

While the extended Higgs sector does alter the results 
from the single Higgs double scenarios, the broad fea- 
tures of the mass splitting structures and preference for 
low masses, particularly for ^4, remain. These features 
are largely a result of the structure of the contributions to 
the electroweak oblique parameters from the fourth gen- 
eration of sequential fermions. Similar results might be 
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expected in almost any extension to the Higgs sector that 
is broadly consistent with a Standard Model-like Higgs. 

IV. SUMMARY 

While stringent limits on m{t^) and 771(64) have been 
found in specific decay modes by the LHC, completely 
ruling out the fourth generation hypothesis requires an 
analysis [43j that combines the results from a number of 
modes to obtain a result that is independent of quark 
mixing in the fourth generation. 

We have used sampling methods to determine the 
probability densities of the masses of a possible fourth 
sequential generation of fermions in scenarios with one 
or two Higgs doublets. With a single Higgs doublet and 
a low (115 or 124.5 GeV) Higgs mass, fourth generation 
mass splitting in the quark sector is less than M\v (sec 
also [11[). Quark sector mass splittings less than M\y 
are favored but less certain if the Higgs mass is 600 GeV. 
A fourth generation is on the verge of being ruled out in 
the case of a single Higgs doublet [3] , but a Type II two 



Higgs doublet model can be designed to reproduce the 
hints at 124.5 GeV from the LHC and the Tevatron. In 
that case, quark mass splittings less than Mw are still 
favored. In all of our scenarios, the most favored values 
for miji) are just above the experimental limit of 110.8 
GeV, making searches for a fourth generation charged 
lepton an interesting possibility. 
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